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The expression and functions of microRNAs (miRNAs) in chronic temporal lobe epilepsy (TLE), the most 
common type of refractory epilepsy in adults, are poorly understood currently. In this study, status 
epilepticus evoked by amygdala stimulation was used to establish rat chronic TLE model. Two months later, 
high-throughput sequencing was employed to investigate miRNA expression profile in rat hippocampus, 
and six miRNAs were confirmed to be differentially expressed. Kyoto Encyclopedia of Genes and Genomes 
pathway analysis indicated that most of the target genes for these six miRNAs were associated with neuronal 
apoptosis. Meanwhile, the levels of miR-423-3p and miR-296-5p were correlated with the activity of 
caspase-3, an apoptosis indicator. Additionally, the loading of miR-423-3p was increased in RNA-induced 
silencing complex whilst caspase-6, a target of miR-423-3p, was reduced in chronic TLE rats. Collectively, 
our findings suggest that miRNAs may exert anti-apoptotic effects in chronic TLE. 

Chronic temporal lobe epilepsy (TLE) is the most common type of partial epilepsy in adults 1 . Symptoms in 
chronic TLE consist of partial or generalized seizures that originate from the amygdala many years after an 
initial brain insult such as status epilepticus (SE), encephalitis or febrile convulsions 2,3 . To date, the 
underlying molecular mechanisms of chronic TLE are still unclear. MicroRNAs (miRNA) are small non-coding 
RNAs that regulate the expression of target genes by binding to the 3' -untranslated regions of target mRNA 4 . In 
mammalians, more than 50% of all miRNAs are expressed in the brain 5 , suggesting a particularly significant role 
for miRNAs in brain physiology 6,7 . In recent years, accumulating evidence suggested that miRNA was implicated 
in the pathogenesis of epilepsy 8,9 . Several studies indicated a differential expression of several miRNAs in brain of 
SE patients and animal models 1012 . Meanwhile, by employing rat model of SE, our group and Hu et al. have 
demonstrated that several differentially expressed miRNAs were involved in modulation of neuronal apoptosis, a 
main pathological characteristics of SE 13,14 . 

Considering the above evidence, we hypothesized that refractory epilepsy such as chronic TLE also led to 
dramatic alterations in brain miRNA levels, and dysregulation of miRNA might be associated with neuronal 
apoptosis and brain damages in chronic TLE. In order to test this hypothesis, we employed high-throughput 
sequencing to investigate the miRNA expression profiles in hippocampus of chronic TLE rats. Meanwhile, Nissl 
staining was used to detected neuronal losses whilst terminal deoxynucleotidyl transferase-mediated dUTP end- 
labeling (TUNEL) and the activity of caspase-3, an apoptosis indicator, were used to assess neuronal apoptosis in 
hippocampus. In addition, the association between differentially expressed miRNAs and caspase-3 activity was 
evaluated, and the loading of miRNAs in RNA-induced silencing complex (RISC) as well as the protein levels of 
their predict targets were also determined here. It should be noted that we induced chronic TLE in rats by using 
amygdala stimulation in this study, as this model could mimic behavioral, electroencephalographic, and mor- 
phological phenomena of human chronic TLE better than chemoconvulsants-induced TLE model 1517 . 

Results 

We obtained 8,839,947 and 9,080,314 raw reads from control group and chronic TLE group, respectively. As 
demonstrated by Supplementary Fig. SI, the most abundant size class was 22 nt, which accounted for 60.97% 
(control group) and 48% (chronic TLE group) of the total reads, respectively. After discarding the sequences 
shorter than 18 nt and low-quality sequences, 8,736,906 and 8,930,938 clean reads of 18-30 nt from the control 
group and chronic TLE group, respectively, were obtained for subsequent analysis. 
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Figure 1 | Confirmation of six differentially expressed miRNAs using 
qRT-PCR at 2 months after amygdala stimulation. The expression of all 
genes is normalized to the levels of U6 snoRNA; Data are presented as the 
mean ± SD; *P < 0.05 versus control group; n = 6 per group. 

The expression patterns of known miRNA in rat hippocampus at 2 
months after amygdala stimulation were indicated by Supplemen- 
tary Fig. S2. We identified novel miRNAs at 2 months after amygdala 
stimulation by using miRDeep software, which was designed to 
uncover novel miRNAs in deep sequencing data as previously 
described 1819 . The algorithm of miRDeep took into account several 
factors to score the probability that a candidate small RNA repre- 
sented a bona fide miRNA: (1) the presence in the deep sequencing 
data of reads corresponding to typical products of miRNA biogen- 
esis; (2) the stability of the putative premiRNA hairpin; and (3) 
homology to previously identified miRNAs. After discarding reads 
corresponding to previously annotated regions, miRDeep analysis 
revealed the existence of a total of 63 novel miRNAs (Fig. S3). 

After comparing the miRNA expression profile in control group 
and chronic TLE group, we detected 25 miRNAs that were differ- 
entially expressed in rat hippocampus. Among these miRNAs, four- 
teen miRNAs were found to be up-regulated whilst eleven miRNAs 
were down-regulated (Supplementary Table SI). Afterwards, six 
miRNAs (miR-455-3p, miR-345-3p, miR-423-3p, novel miR-54, 
miR-365-5p, and miR-296-5p) were screened out based on the cri- 
teria mentioned in "Methods" section, and qRT-PCR was performed 
to validate the alterations of these miRNA. As shown in Fig. 1, at 2 
months after amygdala stimulation, the expression of miR-455-3p, 
miR-345-3p, miR-423-3p, novel miR-54, and miR-365-5p was sig- 
nificantly increased (P < 0.05) while the levels of miR-296-5p was 
markedly decreased in rat hippocampus (P < 0.05). These findings 
were consistent with results from high-throughput sequencing. 

To predict the function of the six differentially expressed miRNAs 
that confirmed by qRT-PCR, the RNAhybrid and miRanda algo- 
rithms were used to obtain a list of potential target genes. Subse- 
quent Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway 
analysis indicated that most of the predicted target genes such as 
Fadd, Fas, Casp, and Bax were closely relevant to neuronal apoptosis 
(Figure 2). Hence, these six miRNAs may play an important role in 
pathogenesis of chronic TLE by modulation of neuronal apoptosis in 
hippocampus. 

We then employed Nissl staining to investigate the neuronal 
survival rate in rat hippocampus. As shown by Fig. 3A, in rats with 
chronic TLE, a dramatic reduction of the neuronal survival rate was 
observed in CA1 and CA3 regions of hippocampus at 2 months 
after amygdala stimulation. To further investigate whether apopto- 
sis was involved in chronic TLE-induced neuronal losses, we per- 
formed TUNEL assay to evaluated neuronal apoptosis rate in rat 



hippocampus. As indicated by Fig. 3B, in hippocampus of chronic 
TLE rats, a significant increase in neuronal apoptosis rate was noted 
in CA1 and CA3 regions at 2 months after amygdala stimulation. In 
addition, we found that the activity of caspase-3, an indicator of 
apoptosis, was significantly higher in hippocampus of chronic TLE 
rats than that of controls (Fig. 3C). These results suggested that 
apoptosis represented part of the mechanisms underlying the 
chronic TLE-induced neuronal losses in hippocampus. 

To evaluate whether the selected miRNAs were implicated in 
chronic TLE-induced neuronal apoptosis, we assessed the correla- 
tion between these six miRNAs and the caspase-3 activity in rat 
hippocampus at 2 months after amygdala stimulation. Pearson's 
correlation analysis showed a positive correlation between miR- 
423-3p and caspase-3 activity in hippocampus (Fig. 3D, r = 0.874, 
P = 0.023). In addition, a negative correlation was found between 
miR-296-5p and activity of caspase-3 (Fig. 3D, r = —0.881, P = 
0.02). These data indicated that miR-423-3p and miR-296-5p might 
be implicated in the regulation of neuronal apoptosis in chronic TLE. 

Subsequently, we investigated whether miR-423-3p exerted its 
function in brain of chronic TLE rats. We selected miR-423-3p rather 
than miR-296-5p mainly because its alteration is more apparent than 
that of miR-296-5p in chronic CLE (shown in Fig. 1). To determine 
whether up-regulation of miR-423-3p led to a corresponding incre- 
ment in loading of the miRNA into the RISC, we immunoprecipi- 
tated Argonaute (Ago)-2, a crucial component of RISC, from 
hippocampus of controls and chronic TLE rats. It should be noted 
that no significant difference was observed in Ago-2 levels between 
controls and chronic TLE rats (Fig. 4A). As demonstrated by Fig. 4B, 
the levels of miR-423-3p bound to Ago-2 were dramatically 
increased in hippocampus of chronic TLE rats when compared with 
those of controls. Meanwhile, western blot analysis indicated a sig- 
nificant reduction in protein levels of caspase-6, a predicted target of 
miR-423-3p (shown in Fig. 2), in hippocampus of chronic TLE rats in 
comparison to controls (Fig. 4C). 

Discussion 

Here, we provided the first evidence that chronic TLE induced by 
amygdala stimulation led to dramatic changes in the miRNA 
expression profile in rat hippocampus. Five miRNAs including 
miR-455-3p, miR-345-3p, miR-423-3p, novel miR-54, and miR- 
365-5p were significantly increased whist miR-296-5p was remark- 
ably decreased in brain at 2 months after amygdala stimulation. 
These observations were in agreement with a previous study from 
Risbud and colleagues, as they found that miR-423-3p was up-regu- 
lated in a chronic TLE model induced by pilocarpine 20 . In addition, 
our findings were also supported by a recent study from Hu et al., as 
miR-296-5p was revealed to be down-regulated in rats with chronic 
TLE 13 . Interestingly, in contrast to our results, the expression of miR- 
455-3p, miR-345-3p, miR-423-3p, novel miR-54, and miR-365-5p 
remained unchanged in a recent study from Gorter and colleagues 15 . 
Meanwhile, Bot et al. found a moderate but significant reduction in 
miR-345-3p in rats with chronic TLE 16 . The lack of reproducibility in 
the miRNA expression profile can be explained by the diversity in 
TLE models, time-point, and the criteria for miRNA screening. Of 
note, the miRNA expression was profiled using microarrays in stud- 
ies from Gorter et al. and Bot et al., whereas we employed high- 
throughput sequencing to characterized miRNA expression here 1516 . 
The differences in the sensibility and specificity between microarrays 
and high-throughput sequencing may also be responsible for these 
contradictory findings. It should be noted that another study from 
our group indicated that miR-345-3p and miR-365-5p were also 
dramatically increased in rats with status epilepticus, implying that 
miRNAs play a crucial role not only in chronic TLE but also in other 
types of epilepsy 14 . Nevertheless, the expression of miR-345-3p and 
miR-365-5p was increased at 24 hour and started to be recovered at 3 
weeks after onset of SE whilst that was up-regulated at 2 month after 
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Figure 2 | Genes potentially targeted by six differentially expressed miRNAs. Boxes in yellow represent miRNAs, whereas cycles in red represent 
mRNAs. Apoptosis-related mRNAs are indicated by blue cycles. 



induction of chronic TLE, possibly attributing to the differences in 
pathophysiology mechanisms underlying these two types of epilepsy. 

As a form of programmed cell death, apoptosis was implicated in 
chronic TLE-induced neuron losses and brain damage 12 ' 2124 . In the 
current study, we found that chronic TLE dramatically increased the 
number of apoptotic neurons in rat hippocampus. These findings 
were supported by a recent study from Hu et al., which demonstrated 
that apoptosis represents part of the mechanisms underlying the 
chronic TLE-induced neuronal losses 13 . Chronic TLE-induced apop- 
tosis can be triggered by activation of intrinsic pathway components 
including pro-apoptotic proteins belonging to Bcl-2 and caspase 
family. In addition, the extrinsic pathway mediated by TNFR1 and 
Fas is also found to be activated following chronic TLE and triggers 
neuronal apoptosis 25,26 . Interestingly, the subsequent KEGG pathway 
analysis indicated that these six miRNAs were targeted to genes 
involved in both intrinsic and extrinsic pathway of apoptosis, such 
as Fas, Casp6, Casp9, and Bax. Meanwhile, we found that the express- 
ion of miR-423-3p and miR-296-5p were significantly correlated 
with the activity of caspase-3, an indicator of apoptosis, further 
implying that these miRNAs participate in the pathogenesis of 
chronic TLE by modulating neuronal apoptosis in hippocampus. 
These findings were in agreement with a previous study from Hu 
and colleagues, which demonstrated expression of several miRNAs 
was dramatically changed and provided modulation on neuronal 
apoptosis in rats with chronic TLE 13 . More importantly, the loading 
of miR-423-3p was found to be increased in RISC where the miRNA- 
based RNA silencing occurred 8 . In addition, the protein levels of 



caspase-6, a predicted target of miR-423-3p, showed expected reduc- 
tion in hippocampus of chronic TLE rats. These results further 
implied that miR-423-3p might be compensatorily up-regulated in 
response to cellular apoptosis, and exerted anti-apoptotic effects in 
chronic TLE. 

In conclusion, we employed high-throughput sequencing to char- 
acterize miRNA expression profile in hippocampus of rats with 
chronic TLE. Six miRNAs including miR-455-3p, miR-345-3p, 
miR-423-3p, novel miR-54, and miR-365-5p were confirmed to be 
differentially expressed at 2 months after amygdala stimulation. In 
addition, miR-423-3p and miR-296-5p were significantly correlated 
with neuronal apoptosis in hippocampus. Moreover, the loading of 
miR-423-3p was increased in RISC whilst the protein levels of cas- 
pase-6, a predicted target of miR-423-3p, showed expected reduction 
in hippocampus of chronic TLE rats. These findings indicate that 
miRNAs may provide protective effects against neuronal apoptosis 
in chronic TLE, and therapies targeting miRNA may open up new 
avenues for the treatment of refractory epilepsy. 

Methods 

Adult male Sprague — Dawley rats weighing 230-270 g were obtained from 
Experimental Animal Center of Qingdao University. They were allowed free access to 
food and water and were maintained on a 12 h light/dark cycle with a controlled 
temperature and humidity, and given free access to food and water. All experiments 
were performed in strict accordance with National Institute of Health Guide for the 
Care and Use of Laboratory Animals. Animal care and sacrifice were conducted 
according to methods approved by the Qingdao University Animal Experimentation 
Committee (permit No. QUEC-130629). 
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Figure 3 | Neuronal apoptosis in rat hippocampus and its correlation with miRNA levels at 2 months after amygdala stimulation. Neuronal survival 
was assessed by Nissl-staining. (A) Representative photo of Nissl-staining in CA1 region of rat hippocampus at 2 months after amygdala stimulation. 
Neurons with dark violet nucleus and intact morphology were identified as surviving neurons. Scale bars: 50 um. The neuronal survival rate was defined 
as follows: neuronal surviving rate (%) = 100X (count of surviving neurons/total count of neurons). Neuronal apoptosis was determined by TUNEL 
assay. (B) Representative photo of TUNEL in rat hippocampus at 2 months after amygdala stimulation. Photos were converted to black and white to 
obtain a better contrast ratio. Neurons with deep black nuclei were identified as TUNEL-positive cells (indicated by red arrows). Scale bars: 50 |im. The 
neuronal apoptosis rate was defined as follows: neuronal apoptosis rate (%) = 100X (count of apoptotic neurons/total count of neurons). (C) The 
activity of caspase-3 in hippocampus was measured using a colorimetric assay kit (Abeam) according to the manufacturer's instructions. (D) Pearson's 
correlation analysis showed that the levels of miRNA-423-3p and miRNA-296-5p were significantly correlated with caspase-3 activity in hippocampus. 
All data are presented as the mean ± SD; n = 6 per group. 
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Figure 4 | Loading of miRNA-423-3p in RISC and protein levels of its predicted targets. (A) Protein levels of Ago-2, a crucial component of RISC, in 
hippocampus of chronic TLE rats and controls were detected by western blot (cropped gels/blots are used here, and the full-length gels and blots are 
shown in Supplementary Fig. S4A). (B) miRNA-423-3p levels elutedbyimmunoprecipitation of Ago-2 in hippocampus of chronic TLE rats and controls. 
(C) The protein levels of caspase-6, a predicted target of miRNA-423-3p, in hippocampus of chronic TLE rats and controls were assessed by western blot 
(cropped gels/blots are used here, and the full-length gels and blots are shown in Supplementary Fig. S4B). (D) Illustration showing that miR-423-3p 
might be compensatorily up-regulated in response to cellular apoptosis and tried to inhibit apoptosis via down -regulation of proapoptotic proteins. All 
data are presented as the mean ± SD; n — 6 per group. 



Rats were fixed in a stereotactic apparatus (Stoelting, USA) under deep anesthesia 
(10% chloral hydrate, 3.5 mL/kg, i.p.). An electrodes were implanted into the right 
basolateral amygdala as previously described (AP: —3.0 mm; L: —4.8 mm; V: 
— 8.8 mm) 27 ' 28 . The electrode was connected to a miniature receptacle, which was well 
embedded in the skull with screws and dental acrylic cement. After the surgery, the 
rats were treated with penicillin for 3 days to prevent infection and allowed to recover 
from surgery for 10 days. 

Seizure was induced by 20 min amygdala stimulation (100 ms train of 1 ms, 60 Hz 
bipolar pulses, 400 uA, every 0.5 s) using a ML1 101 electronic stimulator. 
Electroencephalograms of the right amygdala were recorded with a digital amplifier 
(AD Instrument, USA). Seizure severity was assessed according to a 5-stage scale 29 : (I) 
facial movement; (II) head nodding; (III) unilateral forelimb clonus; (IV) bilateral 
forelimb clonus and rearing; and (V) bilateral forelimb clonus and rearing and falling. 
Following electrical stimulation, rats were video-monitored for 2 months. The rats 
with chronic TLE were identified by occurrence of frequent seizures (at least 2 times 
IV/V spontaneous seizures per week) from 1 week after electrical stimulation 13 ' 15 . The 
time between the last spontaneous seizure and the animals were killed was 2 h. Rats in 
control group were handled in the same manner while did not receive any electrical 
stimulation. 



At 2 months after amygdala stimulation, rats with chronic TLE were anesthetized 
with 10% chloral hydrate (3.5 mL/kg, i.p.) and were transcardially perfused with ice- 
cold saline. Afterwards, brains were rapidly removed and frozen in liquid nitrogen for 
subsequent analysis. 

Total RNA was extracted by the mirVana™ miRNA Isolation Kit (Ambion, USA) 
according to the manufacturer's instructions. RNA quality and quantity were eval- 
uated using Agilent 2100 Bioanalyzer. Total hippocampal RNA from rats with 
chronic TLE as well as control rats was pooled separately and used for miRNA 
sequencing (n = 6 for each group). The small RNAs were ligated to a 5' and a 3' 
adaptor sequentially and then converted to cDNA by reverse transcription. The two 
generated small cDNA libraries were amplified by PCR with primers complementary 
to the adaptor sequences. Afterwards, the libraries were deep sequenced directly using 
the Illumina Cluster Station and Genome Analyze (Illumina Inc, USA) at Beijing 
Genomics Institute 30 . The deep sequencing data have been deposited in NCBI SRA 
database and are accessible through GEO series accession number GSE52443 (http:// 
www.ncbi.nlm.nih. gov/geo/query/acc.cgi?acc— GSE52443). The differences in the 
quantities of miRNAs were determined by comparing the log 2 -ratio of the chronic 
TLE rats and control rats copies. A miRNA was considered "differentially expressed" 
only if it met the following criteria: (a) having at least 10 copies in either chronic TLE 
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or control group, (b) fold-change (log 2 ) > 1.5 or fold-change (log 2 ) < — 1.5 between 
these two groups (P < 0.05 and P 9^ 0). Of note, the Benjamini and Hochberg multiple 
testing correction was applied to adjust the P-values 31 . 

To validate the initial results of high-throughput sequencing in this study, we 
selected miRNAs that met the following criteria for additional qRT-PCR analysis: (a) 
having at least 30 copies in either chronic TLE or control group, (b) Showing fold- 
change (log 2 ) > 1.5 or fold-change (log 2 ) < — 1.5 between these two groups (P < 0.01 
and P 7^ 0). As a result, six differentially expressed miRNAs were screened out for 
subsequent qRT-PCR confirmation. The detailed protocol of qRT-PCR was described 
as previously 14 . All experiments were performed in triplicate. 

Potential target genes of the five known miRNAs were predicted using the database 
of RNAhybrid and miRanda as previously described 32-34 . Meanwhile, the potential 
target genes of novel miR-54 were predicted using a method described by Barozai 
et al 35 . and Wu et al 36 . First, we submitted the novel miRNA sequences as queries to 
the NCBI Basic local alignment search tool program with the following parameters: 
(1) database reference: RefSeq RNA; (2) organism: Rattus norvegicus; and (3) pro- 
gram selection: highly similar sequences (megablast). Afterwards, the mRNA 
sequences showing 75% query coverage were selected and submitted to RNAhybrid 
for the confirmation of the potential targets. Only targets having stringent seed site 
located at positions 2-7 from the 5' end of the miRNA along or with the supple- 
mentary site were selected and further submitted to the miRanda to reconfirm the 
RNAhybrid results. Additionally, the KEGG pathway database was used to filter the 
enriched pathways of miRNA targets. 

Nissl staining was used to detect surviving neurons in hippocampus. Briefly, the 
frozen slices were stained with 1% cresyl violet at 50 °C for 5 min as previously 
described 37 , and the neurons with dark violet nucleus and intact morphology were 
identified as surviving neurons. Six fields of CA1 and CA3 regions in hippocampus of 
each section were randomly selected, and the number of surviving neurons was 
counted. The neuronal survival rate was defined as follows: neuronal surviving rate 
(%) = 100X (count of surviving neurons/total count of neurons). Meanwhile, 
TUNEL was performed using a modified protocol described by Stahelin and collea- 
gues to assess neuronal apoptosis 38 . TUNEL-positive neurons with condensate nuc- 
lear were identified as apoptotic neurons 39 ' 40 . Six different areas were selected at 
random from the CAl and CA3 regions of each section, and the number of TUNEL- 
positive neurons in these areas was counted. The neuronal apoptosis rate was defined 
as follows: neuronal apoptosis rate {%) — 100X (count of apoptotic neurons/total 
count of neurons). In addition, the activity of caspase-3 in hippocampus was mea- 
sured using a colorimetric assay kit (Abeam) according to the manufacturer's 
instructions. 

The association between miRNAs and Ago-2, a crucial component in RISC, was 
investigated as described by Henshall's group 10 41 . Briefly, the hippocampus of six rats 
from each group were pooled and homogenized in 0.7 ml of ice-cold immunopre- 
cipitation buffer. The homogenate was centrifuged at 16 000 X g for 15 min at 4°C 
and the supernatant was retained. Antibodies against Ago-2 (1 : 50; Cell Signaling 
Technology Inc., USA) was added to 400 mg of lysate in a final volume of 1 ml, mixed 
and incubated overnight. Next, 20 ml of 50% Protein- A/G-agarose beads was added 
for 1 h at 4 C and then centrifuged at 16 000 X g for 15 min. The pellet was washed in 
immunoprecipitation buffer, and miRNA was extracted. Afterwards, qRT-PCR was 
performed to quantify the expression of miRNA. 

Western blot analysis was performed as previously described 42,43 . Briefly, hip- 
pocampus of each rat was homogenized and the total proteins were extracted by RIPA 
lysis buffer. Different samples with an equal amount of protein were separated on 8- 
12% SDS polyacrylamide gels, transferred to nitrocellulose membranes, and blocked 
in 5% BSA powder in 1 X trisbuffered saline with 0.1% Tween20 (1 X TBST) at room 
temperature for 2 h. Membranes were incubated overnight at 4 L 'C with a mouse 
monoclonal antibody against Ago-2 (1 : 1000; Cell Signaling Technology Inc., USA) 
and Caspase-6 (1 : 800; Cell Signaling Technology Inc., USA), then washed with 1 X 
TBST, and incubated with HRP-coupled secondary antibody for 2 h at room tem- 
perature. After washing, protein bands were detected with chemiluminescent HRP 
substrate (Thermo Scientific Inc., USA) for 5 min at room temperature and exposed 
to an X-ray film. The band intensity was analysed using Quantity One software 4.6.2 
(Bio-Rad Laboratories Inc., USA) and normalized to a loading control |3-actin. 

Statistical analysis was conducted using the SPSS 17.0 software. Statistical differ- 
ences were determined by independent sample t-test. Pearson's correlation analysis 
was used to investigate to the correlation between the miRNA levels and caspase-3 
activity. All data are expressed as mean ± SD. P < 0.05 was considered statistically 
significant. 
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